Ordered states of adatoms on graphene. 
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We show that a dilute ensemble of epoxy-bonded adatoms on graphene has a tendency to form a 
spatially correlated state accompanied by a gap in graphene's electron spectrum. This effect emerges 
from the electron-mediated interaction between adatoms with a peculiar 1 /r 3 distance dependence. 
The partial ordering transition is described by a random bond three state Potts model. 
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Graphene (monolayer of graphite) is a truly two- 
dimensional crystal, just one-atom-thick pQ. It is a gap- 
less semiconductor with charge carriers mimicking rela- 
tivistic dynamics of massless Dirac fermions [2J, a pe- 
culiarity dictated by the bonding of carbon atoms into 
a highly symmetric honeycomb lattice. Graphene can 
host various adsorbents, in particular atoms, retaining 
its own structural integrity. Such chemisorbed atoms 
(adatoms) may strongly affect electronic properties of 
graphene j3j 131 [5j [6] introducing symmetry-breaking per- 
turbations into the lattice. The type of symmetry break- 
ing depends on the position of the adatom in the hexag- 
onal unit cell of the crystal. In particular, alkali atoms 
position themselves over the centres of the hexagons [7] . 
Oxygen, nitrogen, boron, or an additional carbon [8] pre- 
fer 'epoxy' bonded positions (e-type) and reside above the 
middle of a carbon-carbon bond. Atomic hydrogen and 
halogens reside in the symmetric on-site position above 
the carbon (s-type) [9]. It has also been noticed that 
a pair of hydrogen atoms on the neighbouring sites of 
graphene lattice forms a stable H-H dimer which acts as 
an e-type adsorbent [TP] . 

Here we predict that an ensemble of e-type adatoms 
(those perturbing C-C bonds) tend to order, mimick- 
ing a superlattice structure, even when graphene cov- 
erage by adsorbents is low. The underlying mecha- 
nism is a long-range electron-mediated interaction be- 
tween adatoms similar to the RKKY exchange between 
localized spins in metals [11]. The effect is peculiar to 
graphene. Unlike metals, charge neutral graphene has 
a point-like Fermi surface positioned in the corners K 
and K' = — K of the hexagonal Brillouin zone - called 
valleys. The electron density of states vanishes at the 
Fermi level. As a result, the Friedel oscillations in charge 
neutral graphene are commensurate with its honeycomb 
lattice and decay as the inverse cube of the distance to 
the adatom [IS]. We show that such an interaction in a 
dilute ensemble of e-type adsorbents may result in their 
partial ordering associated with a superlattice structure 
with the unit cell three times larger than in graphene, 
as illustrated in Fig. [I] We present our results in the 
following order. Starting with a particular tight-binding 




FIG. 1: An ordered ensemble of adatoms (blue) on graphene 
lattice: all impurity atoms occupy positions at the intersec- 
tion of graphene bonds and the bonds of the underlying ficti- 
tious superlattice. The inset shows how the superlattice orig- 
inates from the scattering matrix W, Eq. ([3]). To each bond 
there corresponds a pair of vectors (u, v) forming a periodic 
pattern with three times graphene lattice's period. 



model for an e-type adsorbent, we determine the form of 
a perturbation it creates for the electrons in graphene. 
Using group theory we classify such interactions beyond 
a specific microscopic model and determine the condi- 
tions under which RKKY interaction between adatoms 
leads to a partially ordered state with a gapful electronic 
spectrum. We conclude by discussing experimental sig- 
natures of the effect. 

The 7r-electron band in graphene is well described 
by the closest-neighbor tight-binding model, H = 
7oE[r A r s ]{4 A c r B + h.c.} with hoping parameter 70 ~ 
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TABLE I: Symmetry-based classification of RKKY interactions for different types of adatoms. The parameterization of the 
adatom scattering matrix W a is given in terms of an orbit Wj of the lattice symmetry group in a given irrep i. The potential 
f2„ 6 , is a contribution of the given orbit to the adatom-adatom interaction, Eqs. |5|. The last column indicates ordering favored 
by the given interaction term. 



Irrep i 


Orbit Wi 


ab 


type: position 


Partial order 


Ax 


1 




All 




B 2 


sE z A 2 , s — ±1 






Sublattice 




V ( A 11*1 ii a XT 

^z(A • UJ, U G V3 


Uo • U(, 




3 x unit cell 
superlattice 


E 2 


A z (£-v), vG V 3 


v • vj, — |(n • v a )(n ■ Vb) 




None 


G 


(v x S)(u x A), u,v G V 3 


i(u a • u 6 )[3(n ■ v a )(n • vj,) - (v a ■ v b )] 




None 


E'l 


sE x (A ■ su) + E S (A x su) u G V 3 , 
s = diag(l,s) ' s = ±1 


2 ~ 3 2 aS '' [(n • u a )(n ■ u b ) - s a s b (n x u a )(n x u b )] 




None 



by the three matrices A x . y = cr| <g> <rj and A z = I s ® <rj , 
satisfying [Aj,£j] = 0. This symmetry manifests of the 
conservation of the electron's valley index. Matrices £, A 
and their products AjXj can be arranged into irreducible 
representations (irreps) of the symmetry group Q of the 
honeycomb lattice [T5J [T5] , which includes lattice transla- 
tions, Cq v rotations and mirror refections. All operators 
Hi and A,,- change signs upon time inversion, therefore 
only products AjEj are time-inversion-symmetric and 
can appear in the scattering matrix [T7] W a , Eq. ^ , de- 
scribing static perturbations |T5J [TS] created by adatoms. 

For an adatom of a general symmetry type W a , can 
be expanded into orbits in the irreps [SO] of G, W a = 

\iWi- The classification of orbits by the irrep and the 
symmetrytype of adatom is given in the second column 
of Table |l| In particular, the matrix W a of an e-type 
adatom (or the H-H dimer) is expanded into orbits in 
the irreps Ai,E[,E2 and G, corresponding to the four 
terms in Eq. ^ . Vectors u and v in Eq. (|3| take values 
in the set V3 shown in Fig. [T| three unit vectors on the 
x-y plane, at 120° angles. Each given bond of graphcnc 
lattice is characterized by a pair (u, v) as shown in Fig. [T| 
The distribution of u and v forms a periodic pattern with 
three times the graphene lattice period: the intervalley 
scattering implies the momentum transfer AK = K — K' 
such that 3AK is a reciprocal lattice vector. Thus, W a in 
the e-type case is periodic on a superlattice, whose unit 
cell is three times as big as graphene 's and contains nine 
distinguishable C-C bonds. The three non-trivial terms 
in Eq. ^ resemble the coupling of an electron to the in- 
plane T- and K-point phonons [16] - Indeed, the E[ term 
parameterized by the vector u resembles the effect of the 
K-point breathing phonon mode. Below, the ensemble 
average of u will play the role of the order parameter. 
The G term is similar to the 4-fold degenerate K-point 



3eV 0. The sum runs over all pairs [r^rg] of neighbor- 
ing A and B sites of the lattice and c^/c are the on-site 
electron creation/annihilation operators [T3]. An e-type 
adatom attached to the bond between the cites r' A and 
t' b creates a local perturbation 

H a = 7ofi(c^c r;4 +c^c r ' B )+7o6(c^c r ' B +c^c r ^), (1) 

where £1 , £2 < 1 determine how the adatom affects the 
on-cite potential (£1) and the electron hopping amplitude 
between the cites (£2)- 

The long-range RKKY interaction between two 
adatoms is due to the perturbation of the electron 
spectrum near the Fermi energy and is adequately de- 
scribed in terms of the four-component field, tp(r) — 
(ipi, ip2, ">p3, ipi) , which is smooth on a scale of the lattice 
constant, o: 

= (8)J B Je flCr W') + e" fKr Wr) 1 r = r A 
U) \ e iKr V> 2 (r) + e- iKr V> 3 (r), r = r B ' 

In the presence of an adatom ip(r) obeys the Hamiltonian 

m 

H = v f^(E-p)ijd 2 r + rlva^(r a )W a iP(r a ), (2) 

W a = A Al +A B ;S Z (A-U a ) 

+ A E2 A 2 (E • v„) + A G (A x u a ) • (53 x v a ), (3) 

with v — v / 3a7o/2?i and T, XtV — o% y ® a\. The Pauli 
matrices cr", a y and cr" operate on the valley (a = v) or 
sublattice (a = s) indices. Together with S z = a\ ® l v , 
matrices £ X)3/ form a representation of the SU2 algebra. 
The first term in H determines the Dirac electronic spec- 
trum. It possesses a "flavour" SU^ symmetry generated 
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lattice mode. The Ei term resembles the uniaxial strain 
due to a T-point optical phonon. The parameters Xi in 
W a are specific for particular atoms. For the model in 
Eq. ([I), X E[ = X E2 = £ 2 and X Al = X G = 

Unlike e-type adatoms, the perturbation introduced by 
an s-type adatom (e.g., H on a lattice site) cannot be 
related to the in-plane phonons, since it explicitly dis- 
tinguishes A and B sublattices of graphene. The two 
orbits encountered in the matrix W a in this case can be 
related to the out-of-plane phonons in the presence of a 
transverse electric field (z — > — z asymmetry): B% resem- 
bling the T-point and E'{ the K-point phonon. In both 
cases, the A/B-residency of an adatom is accounted for 
by s = ±1. 

Each type of the adsorbents listed in Table |T] creates 
Friedel oscillations of the electron density breaking the 
symmetry in the same way as the adatom does. The po- 
larization caused by one adatom extends over long dis- 
tances, thus leading to interaction between the adatoms. 
The energy of this interaction between a pair of adatoms 
at a distance |r a b| ^> a can be expressed through the 
imaginary time (Matsubara) Green function of electrons 
in a clean graphene, G(r, r) [2T] as 
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Here, trace is taken over the valley and sublatticc indices, 
and spin degeneracy is taken into account. Equation Q 
yields a long-range pair correlation energy, 
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Interaction between adatoms of the same type depends 
on both the type of the adatoms and the parameters 
Xi characterizing the adatom-electron coupling in each 
symmetry- breaking interaction channel Wf. 



^ab - ^ 



(6) 



For the model in Eq. ^ with Xe[ larger than other cou- 
pling parameters, fl ab « X 2 E ,u a ■ u b , so that the inter- 
action between lateral degrees of freedom of the two ad- 
sorbents looks like an isotropic ferromagnetic exchange. 
Therefore the adatoms tend to occupy preferably bonds 
with the same u and form a partially ordered state shown 
in Fig. [l] The transition to such a state is described by a 
special case of random exchange three-state Potts model 
[22] : 'spins' u reside on randomly distributed with den- 
sity p sites and experience pairwise exchange interaction 
— Ju a -u b /r 3 . According to recent cluster Monte-Carlo [23 
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FIG. 2: The schematic phase diagram of the epoxy-bonded 
adatom system on graphene. The critical surface separates 
the gapped ordered phase (Fig.[T]) from the disordered phase. 



studies this model undergoes an order-disorder transition 
at a critical temperature T c ss 8p 3 / 2 J. For the Eq. ([lj 
model T r is evaluated as 



0.6X 2 E ,(a 2 P f/ 2 



hv 

a 



(7) 



Consider now effects of other terms in ^ allowed by 
the e-type symmetry. These terms are listed in the rows 
Ai , E 2 and G in Table |T] The symmetric perturbation 
parameterized by X Al leads to the repulsion between 
adatoms regardless of which bonds of the extended su- 
percell they occupy. The coupling parameterized by Xe 2 
causes an anisotropic 'antiferromagnetic' interaction of 
the alternative set of Potts "spins" v a . Frustration pre- 
cludes ordering of v Q , which could lead to a unilateral de- 
formation of the lattice. However, a presence of X% 2 fl^ 2 ^ 
in fl ab , Eq. ^ does not affect the ordering of vectors u a , 
at least up to the quadratic order in Xe 2 - The isotropic 
interaction leading to the ordering only competes with 
the anisotropic interaction between adatoms caused b 
the last term in W a , Eq. ^ (the fifth row of Table |l[ 
If the latter is strong, Xq > A^j, it suppresses order- 
ing by frustration. As a result, T c of the order-disorder 
transition decreases with increasing ratio Xq / X E ^ until it 
vanishes at a quantum critical point. 

Increasing the density n e of mobile carriers in graphene 
should also suppress T c . Indeed, at finite n e the RKKY 
interaction develops Friedel oscillations [19], F oc (u a ■ 
Ub) sin 2 (fcp r)/r 2 , which lead to a random sign of the 
exchange coupling between adatoms at distance r ab > 
y/l/n e . At sufficiently large ratio n e /p this effect should 
completely destroy ordering of Potts "spins". These 
considerations are illustrated by the phase diagram in 
Fig. [2] where the "quantum critical line" corresponds to 
the parametric condition T c (n e /p, Xq/X E ') — > 0. Further 
analysis of these transitions is beyond the scope of this 
paper. 
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The interaction of s-type adsorbents residing on the 
honeycomb lattice sites (e.g., hydrogen) is described in 
rows 2 and 6 of Table [i] Each s-type adatom can be char- 
acterized by the Ising "spin" s a taking the value +1 or 
— 1 depending on which sublattice A or B it occupies. 
These spins may, potentially, establish sublattice order- 
ing. However, there will be no ordering of the "spins" 
u a : the interaction described in Table [I] is anisotropic 
and causes frustrations. In contrast, a pair of hydro- 
gens forming an H-H dimer on the nearest A/B sites |10j 
falls into the same symmetry class as e-type adatoms and 
can establish the same type of ordering. It has been no- 
ticed that in hydrogenated graphite both configurations 
of H atoms are present [10] . Since the decomposition of 
the interaction in Eq. ([5| suggests the absence of mu- 
tual correlations between adsorbents of the e-type and s- 
type, we expect ordering of the H-H dimers on graphcnc, 
even if only a fraction of hydrogens covering the flake is 
dimerised. 

The predicted ordering will strongly influence trans- 
port and optical properties of the material. As tempera- 
ture approaches T c from above large clusters of ordered 
phase will act as intervalley-scattering Bragg mirrors for 
electrons. Back-scattering from such mirrors should lead 
to a power-law increase of resistivity near T c |25j . At 
T < T c , the spectrum of electrons becomes gapful, 

s = ±^v 2 p 2 + A 2 , A= — A^pa 2 »T c , (8) 

a 1 

as can be seen from the mean-field Hamiltonian |24j . 
(H) = vll ■ p + hvpa\ E 'Y, z K ■ u, valid for p < hp 1 / 2 . At 
low carrier density this will lead to the activated trans- 
port regime typical of semiconductors. 

Partial ordering of e-type adsorbents should also be 
manifest in the structure of the D-peak in the Raman 
spectrum. The D-peak is associated with the excita- 
tion of one optical phonon at momentum K (or K') 
and is forbidden by momentum conservation in pristine 
graphene. In the presence of random scatterers, the D- 
peak is seen as a low-intensity, I(T > T c ) cx p\ 2 E , , fea- 
ture strongly broadened due to the disorder-induced un- 
certainty in the emitted phonon momentum. Domains 
of ordered adsorbent, with the size L, will scatter elec- 
trons between valleys K and K' coherently. This will 
enhance the intensity of the otherwise forbidden transi- 
tion, I(T < T c ) oc p\%, x pL 2 > g(T > T c ), and re- 
strict the uncertainty of the emitted phonon momentum 
to 5q l/L. Thus, one can predict that the ordering 
of adatoms abruptly enhances and narrows the D-peak 
in the Raman spectrum. To mention, an observation of 
a hopping conductivity accompanied by a sharp high- 
intensity D-linc in the Raman spectrum in graphene ex- 
posed for a long time to hydrogen atmosphere has been 
reported in Ref . . We predict a similar behavior of the 
ARPES spectrum of graphene: ordering should strongly 
enhance the photoemission of electrons from the center of 



the Brillouin zone at energies close to the Fermi energy. 
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